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We generate high-dimensional time-bin entanglement using a mode-locked laser and analyze it 
with a 2-photon Fabry-Perot interferometer. The dimension of the entangled state is limited only 
by the phase coherence between subsequent pulses and is practically infinite. In our experiment a 
pico-second mode-locked laser at 532 nm pumps a non-linear potassium niobate crystal to produce 
photon pairs by spontaneous parametric down-conversion (SPDC) at 810 and 1550 nm. 



I. INTRODUCTION 



Entanglement is one of the most useful resources for quantum information Q . Most entanglement based experiments 
involved 2-level or eventually 3- level systems However over the last few years systems with higher dimensions 

[JL^ , have received increasing attention for a variety of reasons. The tolerance to noise of quantum key distribution can be 
increased thanks to high-dimensional systems High-dimensional entanglement allows for the required efficiency 
| of detectors to close the detection loophole in EPR experiments to be reduced [||. Moreover, their properties differ 
from classical ones more than 2-levels system and they have greater robustness against noise 0, @ • 

High-dimensional systems can be obtained in two ways. Firstly, we can get multi-photon (more than 2) entanglement 
by using high-order parametric down-conversion JtJ llfj. ITlj . Secondly, we can consider two-photon entanglement in 
high-dimensional systems. This second approach has the experimental advantage of higher coincidence count rates as 
you have to create and detect only two photons. For example entanglement of higher order angular momentum states 
1 of photons has been demonstrated [l2lll8T |. However, time-bins seem to be the ideal scheme for higher dimensional en- 
tanglement. Indeed, mode-locked laser can easily produce entangled states of almost arbitrarily high dimensions. This 
I/"} i has been shown using Michelson interferometer |14L Il5| . i.e two-dimensional analyzers. Unfortunately, the extension 
of this analysis to higher dimensions, using e.g. interferometers with n different paths, dramatically complicates the 
experimental task. In this paper, we present an experimental realization of a high-dimensional analysis using Fabry- 
O ,' Perot like interferometers. We start with a theoretical description of our analyzer before presenting the experiment 
and the results. 
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A mode- locked laser is used to pump a non-linear crystal in order to produce time-bin entangled photon pairs. We 
consider a Z?-pulse train and assume a pair creation probability much lower than per pulse to reduce the creation 
of two pairs in a Z?-pulse train. When a photon pair is created in time-bin j, the state is \j, j). As the time-bin in 



which the photon pair is created is uncertain, the state after the non-linear crystal is of the form: 

D 



^crystal) = ^ CjB^ \j,j) (1) 
3=1 

where Cj are the probability amplitudes and tf>j are the phase difference between successive pulses. For a mode-locked 
pump laser Cj and <pj are constant. 



To analyze the high-dimensional time-bin entangled state we use Fabry-Perot like interferometers (see Fig. QJ. 
During one turn through the interferometer a photon is delayed exactly by one time-bin At = — - — where fi aS er is 
the repetition frequency of the laser. Let us first consider the detectors D a and D . After the interferometers where 
the two photons go to detectors D a and D b , respectively, the state \j,j) evolves as follows (the first passage through 
the interferometer to the detectors adds only a global phase, which is not taken into account): 
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FIG. 1: Simplified scheme of the experiment. Pulses of a mode-locked laser are sent through a KNbOs non-linear crystal 
(NLC) to produce photon pairs. The two output modes are coupled into separated fibers and go through interferometers with 
optical paths corresponding to the distance between two successive pulses. 
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where <p x is the phase applied on the photons in mode x — a, b and can be considered constant for successive turns. 
The first row represents the situation when photon a covers one more turn than photon b, second row represents 
photons a and b in the same time-bin and the third row represents the case with photon b covering one more turn 
than photon a. 

The state of Eq. (JTJ then evolves, according to J3J), to: 

D-l D-n 

\tp) = EE C n . D (e in ^ + e ii{n+l)4>a+<t> b ) + . . . + e i(D<j, a + (D-n)<j> b )^j + 
n=0 j=l 
D-l D-m 

+ E E c m ,D(e vm<t, » + e i(^+( m - 1 )« + ... + e l U D ~ m ^+ D ^)\j,i + m) (3) 

m— 1 j — 1 

To analyze the system we measure the difference in the time of arrival of photons a and b at detectors D a and 
respectively. For each value of n or m from Eq. @ there is a corresponding peak in the histogram of the time 
difference At = tb — t a . The central and highest peak corresponds to coincidences with A< = 0, while the first peak 
on its right corresponds to At = At and first peak on its left to At = — At (see Fig. 01. 



The relative height of these peaks, i.e. the probability of coincidences between D a and Di, for different At, can be 
calculated as follows (for D — ► oo and without losses): 



Pn=0 = Po ~ (tla,tlbt2at2b) 

P n <o = (r 2a ria) 2ln] Po 
Pn>o = (r 2b r lb ) 2n P Q 



where: 



1 



1 - r2ar2bri Q n b e i W<*+« 



(4) 



• P n is the coincidences probability between D a and Db for the n-th peak in the arrival time difference histogram. 
By convention we denote the peak corresponding to photons doing the same number of turn in each interferometer 
by n = 0, the first peak to the right is denoted by n — 1, and so on and similarly with peaks on the left n = — 1 
for the first one and so on. 
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FIG. 2: Coincidences as a function of the difference in arrival time for the two photons at detectors D a and Dt- These 
correspond to the sum of different interfering terms which are due to the different possibilities for each possible arrival time 
difference. 



• tmx (i"mx) is the transmission (reflection) amplitude of the coupler m = 1 (2) for the first (second) one on the 
way of the photon in interferometer in mode x = a, b. By convention a "reflected" photon stays in the same 
fiber. One has to choose t mx -C r mx in order to have a strong weighting of the terms that involve many turns 
in the interferometers, which are characteristic of high-dimensional entanglement. 

For all P n terms the phase dependence is the same for all n so the different peaks of coincidences in the gate of 
detection on D b show synchronous oscillations as a function of the sum of the phases in the interferometers, <f> a + <f>b- 
It's interesting to also calculate the probability to have coincidences between the detectors D a and the third detector 
D' b that we use as a control (see Fig.^l: 
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where: 



• P' n is the coincidence probability between D a and D' b for the n-th peak in the histogram of arrival time difference. 
As for P n we use the convention that n = for the case when photons in modes a and b go to the detectors 
with the same number of complete turns, n = — 1 for the first peak on its left and so on and n = 1 for the first 
peak on its right and so on. Note that this histogram is asymmetrical, as all peaks for n > are much smaller 
than those for n < 0, since ti b ri b . 

The P' terms show the same behavior as the P terms, however we have a minimum of coincidences with detector 
D' b when we have a maximum of coincidences with detector D b as can be expected by conservation of energy. Indeed, 
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FIG. 3: Simulation of normalized coincidences as a function of the phase. The solid line corresponds to coincidences between 
D a and Df, and the dashed line corresponds to coincidences between D a and D' b (D — > oo, r mx — V0.9). 

the light has to go out of the interferometer b by one of the two outputs if there is no losses (absorption) in the 
interferometers. This different behavior for the two terms is due to the minus sign in front of r\ h in the formula of Pq 
which is a consequence of the ^ phase acquired when a photon is "transmitted", i.e. coupled. In Fig. [3] we can easily 
verify that these probabilities sum up to unity in the case without losses. We see that as a function of <p a + <Pb, we do 
not obtain a sinusoidal variation as we are used to in the case of qubits. The curves remind us of the transmission 
through a Fabry-Perot interferometer, which is a consequence of the high dimensionality of interferences, i.e. many 
interfering paths. The goal is now to find this signature experimentally. 



III. EXPERIMENT 

A mode-locked, frequency-doubled Nd-laser (Time-Bandwidth GE-100, A = 532 nm, FWHM< 10 ps, Fi aser = 430 
Mhz, Pmean = 30 mW) is the heart of our experiment (see Fig.0J. A f% = 200 mm achromatic doublet lens focalizes 
the light on a potassium niobate non- linear crystal (KNbOs, Castech, 9 = 23°, ip — 0°) cut in order to obtain collinear 
signal and idler at 810 nm and 1550 nm wavelengths by type I parametric down-conversion. A dichroic mirror is used 
to separate the two non-degenerated photons. In each output arm, a lens is firstly used to collimate the beam while 
the second one focuses light into the monomode optical fiber at 810 nm and 1550 nm respectively. As usual we also 
have to be very careful to filter out all photons originating from the pump. First we remove the remaining photons at 
1064 nm, using a KG5 filter, a dispersive equilateral prism and a diaphragm. In order to remove the pump photons 
at 532 nm after the crystal, we put a RG-610 filter coated with a dielectric mirror at 532 nm and a 10 nm (FWHM) 
bandpass filter at 810 nm in arm a and a combination of an AR coated silicon filter and a 20 nm (FWHM) bandpass 
filter at 1550 nm in arm b. 

The interferometer b is made of two R/T = 90/10 couplers which are spliced together to the required length. An in-line 
fiber polarization controller (Newport PolaRite F-POL-IL) is added in this loop. The realization of interferometer a is 
different to simplify alignment with the interferometer b (see Fig.0J. It is made from a monomode fiber at 810 nm of 
about 23.6 cm length with dielectric mirrors deposited on the cleaved extremities with reflectivity and transmitivity 
R/T = 90/10. The fiber is cut slightly shorter than it normally should be and then it is stretched with a translation 
stage. A piezoelectric actuator (PZT) allows us to then vary the length, i.e the phase, by a few wavelengths. Both 
interferometers are enclosed in separated PI (proportional and integral parameters) temperature-regulated boxes. 

Alignment of the interferometers is the first experimental problem. The optical path lengths of interferometers 
a and b must be the same to within the coherence length of the photon pairs, i.e 120 /im, as well as the cavity 
length of the pump laser, to within the coherence length of the pump (~ 2 mm in fiber). We use an auxiliary, bulk 
Michelson interferometer where the path length difference is firstly adjusted to interferometer b, using low coherence 
interferometry. We then adjust the cavity length of the laser and of the interferometer a to the auxiliary Michelson 
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FIG. 4: Experimental scheme. The pulses of the mode-locked laser are sent through a KG5 filter, an equilateral prism (EP) 
and a pinhole (P) to be monochromatic. After that they go through a KNbOs non-linear crystal (NLC) and produce non- 
degenerate photon pairs at 810 and 1550 nm by type I parametric down-conversion. The photons at 810 nm are transmitted 
through a dichroic mirror (DM) and the photons at 1550 nm are reflected. The photons in the 810 nm arm are filtered by a 
reflector at 532 nm deposited on a RG-610 filter (R + RG filters) and a bandpass filter of 10 nm (FWHM) centered at 810 nm 
(BP a ). The photons in the 1550 nm arm are filtered by a Si and a bandpass filters of 20 nm (FWHM) centered at 1550 nm 
(BPt,). The photons are coupled into monomode (at their wavelengths) fibers and go through the interferometers. Detection 
on the Si APD D a triggers (not represented) the InGaAs/InP APD' s (Df, and D' b ) and starts the TDC. Detections on the 
InGaAs/InP detectors stop the TDC. 



interferometer. 

The 810 nm photons are detected by a silicon (Si) single photon detector in passive mode (EG&G PQ-F830) and the 
1550 nm photons are detected by a InGaAs/InP single photon detectors (ID Quantique, id 200 SPDM) gated by 
the Si detector. The gate width is 50 ns so we can detect 20 time-bins in each gate. The Si detector output starts 
the Time- To-Digital Converter (TDC, ACAM AM-F1) and one stop is given by each output of the two InGaAs/ InP 
detectors Dt and D' b . 



IV. RESULTS 



Figs. 0(a) and (b) show typical histograms for the time difference between a click from detector D a -Db and D a -D' b , 
respectively, recorded with the TDC. In the case (a) the number of accumulated coincidences is much lower than in the 
case (b) , because most of the light is reflected on the first coupler of interferometer b and goes directly to the detector 
D' b . The vertical lines represent the different time windows used in the measurements of Fig. 0. In all measurements 
we record the number of coincidences as a function of <fi a , which is a function of the PZT voltage. For this purpose, 
we accumulate the number of coincidences, typically for 1 minute, then increase, step by step, the voltage on the 
PZT. In order to minimize fluctuations due to varying pump power or coupling of the down-converted photons into 
the fibers, we normalize all coincidence rates with respect to the average single count rates of detector D a . We also 
subtract the noise of InGaAs/ InP detectors, whereas the dark-count of the Si detector (D a ) can be neglected. The 
noise of the InGaAs/ InP detectors is due to the thermal dark count of the detector and it is of 15.6±0.1 Hz on 
and 17.6±0.1 Hz on D' b for an efficiency of detection of about 16 and 18% respectivly, for the entire gates of 50 ns 
and a gating frequency of about 4.6 kHz. 
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FIG. 5: Net coincidences as a function of the difference of arrival times. The histograms are the accumulation of a long term 
measurement, (a) Histogram between detectors D a and D^. The vertical lines represents the different coincidence gates for 
Fig- El f° r the complete gate (solid line), the three central peaks (dashed line) and the central peak (dotted line), (b) Histogram 
between detectors D a and D' b . In the case (a) the accumulated coincidences is lower than in case (b) because T R. 



This gating frequency correspond to the rate of detection on the Si detector and it can be explained as follows. 
The repetition frequency of the laser is 430 MHz and the incident power on the crystal is approximately 17 mW and 
we have a probability of pair creation lower then 1%, thus we have a pair creation frequency of <4.3 MHz. We can 
expect a global coupling factor of the order 10 % between the crystal and the monomode fiber at 810 nm (including 
the losses through the filters). Therefore 430 kHz of photon at 810 nm are coupled into the fiber. With the losses of 
about -14 dB when the light goes through the interferometer a and a detection efficiency for the Si detector of the 
order of 40-50% we can expect a detection frequency of few kHz. 



Firstly, we observe that the coincidences between D a and Db vary synchronously for all different detection windows. 
In Fig.[S]we see the coincidences accumulated over the entire gate of 50 ns, the three central peaks and only the central 
peak, respectively. These three different coincidences sets oscillate synchronously as expected. 

We notice that the peaks are considerably broader than what we would expect for the ideal case according to Eq. Q 
and depicted in Fig. |3| Of course the experiment is not perfect and we can improve our theoretical model in order to 
take into account the following four experimental limitations: a) there are losses in the order of 5% per round trip for 



7 




Time [min] 

FIG. 6: Net and normalized coincidences between Si detector D a and InGaAs / InP detector D), as a function of time while 
changing the phase (see text). We see interferences for the central peak (■), the three central peaks (A) and the entire gate 
(about 20 peaks, ♦) (see Fig.^J. The lines are only presented as a guide. 




FIG. 7: Coincidences as a function of time (see text) with two different fixation systems for the fiber in interferometer a for 
the entire 50 ns InGaAs/ InP detectors gate. In the first case the fiber is squeezed within a metallic holder (■) and in the 
second case the fiber is glued on a metallic holder (♦). The first measure is renormalized to have the same amplitude and time 
dependance. 



interferometers a and b. The losses essentially reduce the contribution for the cases where both photons make several 
roundtrips in the interferometers, b) High visibility interferences can only be achieved if the polarization states of the 
interfering paths are identical. For this purpose, we inject in the interferometer, light from an external and pulsed 
laser, polarized in the same direction as the down-converted photons. We introduce then an auxiliary polarizer at the 
output of the interferometer. With the internal polarization controller, we now maximize the transmission of all the 
peaks, corresponding to zero, one, two and more roundtrips in the interferometer. Unfortunately, a perfect alignment 
is very difficult to achieve in practice and the remaining misalignment increases with the number of roundtrips. In the 
interferometer a we don't have a polarization controller, we count on the fact that in a short straight fiber, without 
stress induced birefringence, the polarization is not altered, in principle. However, we have to pay attention to stress 
induced birefringence. In particular we realized that the peaks are narrower if we glue the fiber on a holder rather 
than fixing it by squeezing it in metallic holders (see Fig. 0) . c) We have to take into account that our light is not 
monochromatic and hence the phase is not exactly the same for all wavelengths, d) Moreover, a slight fluctuation 
in temperature during the measurement can introduce some phase noise in the order of tt/8 per 0.01°C. Again, this 
phase noise adds up with each roundtrip and is hence more important for the terms characterizing the high-order 
entanglement. 

We try to take into account these experimental limitations. We consider the losses as mentioned above and also the 
effective spectrums. The most limiting bandpass filter is the 20 nm at 1550 nm one and it correspond to a bandpass 
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FIG. 8: Measurements and simulations of coincidences as a function of time (see text) for the entire 50 ns InGaAs / InP 
detectors gate. The measured coincidences between detector D a -Df, (D a -D' b ) are represented by ♦ (■). For the simulations we 
consider spectrums of 5.4 nm at 810 nm and 20 nm at 1550 nm and gaussian fluctuations on the phase of about n/8 FWHM 
corresponding to a length fluctuations of ±25 nm of interferometer a. 



filter of 5.4 nm at 810 nm. Finally we introduce gaussian phase fluctuations of tt/8 (FWHM). We can qualitatively 
reproduce the measured curves (see Fig. 0. Hence, we conclude that we have demonstrated the generation and the 
detection of high order entanglement. 



V. CONCLUSION 



We experimentally demonstrated high-order time-bin entanglement. Whereas the creation of high-dimensional 
states is very convenient with mode-locked lasers, we have to realize that the experimental difficulties for the detection 
increase significantly with the dimension of the Hilbert space. Therefore, despite some potential advantages of high- 
dimensional entangled states discussed in the introduction, these states tend to be of limited value for near future 
practical applications. 
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